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abstract 

Multistage cancer nslc models that account for clonal growth of intermediate, 
premalignant cell populations have been successfully used to describe age-specific 
incidence trends for human and experimentally induced animal tumors. Currently there is 
keen and growing interest in applying "cell-kinetic multistage" (CKM) models, which 
distinguish between mutations and cell population kinetics as separate processes 
influencing the rate of tumor formarion, to the problem of predicting human cancer risk that 
I i may be posed by environmental agents. This paper illustrates the application of such 

1 I models to cancer risk prediction for carbon tetrachloride (CCl,)), taking into account rodent 

' cancer-bioassay data and CCU pharmacokinetics in rodents and humans. The increased 

1 importance of knowledge concerning exposure scenarios when using CKM models for risk 

extrapolation, as opposed to more traditional approaches, is illustrated. The presumption, 

' made in most CKM models considered to date, that preneoplasiic cell growth is exponential 
over time, is critically evaluated along with the impact on predicted CCl 4 -induc^ cancer 
risk of using a more biologically realistic cell-growth model. Because the potential 
application of CKM models to cancer risk prediction could have important consequences on 
j public health measures directed at controlling environmental exposure to chemicals, the 
biological plausibility of specific CKM models considered for application deserves renewed 
scrutiny. 

KEYWORDS; Carbon tetrachloride (Ca 4 ), cell-kinetic multistage (CKM) models, phar¬ 
macokinetic (PBPK) models, cancer-risk extrapolation 

i 

I INTRODUCTION 

Multistage canccr-risk models that account for clonal growth of intermediate 
(initiated, premalignant) cell populations, or "cell-kinetic multistage" (CKM) models, have 
long been successfully used to ckscribe age-s^ific incidence trends for human cancer and 
experimentally induced tumors in animals.More recently, a 2-stage stochastic model 
allowing for birth and death of initial and intermediate cells has been used to model human- 
cancer incidence in particular and has been pnroscd to be capable of describing 
spontaneous and induced carcinogenesis in general.^ Currently there is keen and growing 
interest in the possibility of applying this 2-stage model to the problem of predicting the 
human cancer risk posed by environmental agents.^'^^ Because CKM models distinguish 
between mutations and cell-population kinetics as separate processes that influence the rate 
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of tumor formation, such models can predict widely different cancer risks for agents 
affecting only one process (e.g., promotion via induced cell proliferation, assuming a dose- 
response threshold) but not the other (e.g., mutagenesis, assuming a no-threshold dose- 
response relationship). Such an approach allows for distinct treatment of genotoxic 
carcinogens and nongenotoxic tumor promoters, provided the latter can be operationally 
defined. This approach has been advocated^but has not yet been incorporated into 
wddely used cancer-risk-assessment procedures. 

This paper illustrates issues involved in the application of CKM models, in 
conjunction with physiologically based pharmacokinetic (PBPK) models, to the problem of 
cancer-risk extrapolation for the Case of CCI 4 using the following approach. Appropriate 
PBPK models for predicting effective dose in rodents and humans as a function of CCI 4 
exposure are first identified. It is presumed here that CCI 4 is not effective in producing 
mutations assi med (by the CKM approach) to be necessarj' to the carcinogenic process, but 
that CCI 4 is an effective "promoter " by virtue of an ability, at cytotoxic doses, to selectively 
increase the rate of net cell growth within premaJignant ceil clones. Hence, an appropriate 
model is fit to rodent data on the dose-response relationship for CCla-induced cell 
proliferation. This dose-response relationship is assumed to be the same for humans. An 
increase in cell proliferation is then presumed to correspond to a proportional increase in 
the value of the growth parameter of a 2-stage CKM model fit to a rodent-bioassay data set 
on increased cancer risk above background associated with different effective CCI 4 doses, 
which in turn correspond (using the rodent PBPK model) to the exposure scenario used in 
the bioassay selected as the basis for risk extrapolation. 

Increased risk to humans associated with a given effective CCI 4 dose is then either 
presumed to be identical to predicted increased risk to rodents or calculated by fitting a 
corresponding 2-stage CKM mode) to human incidence data for cancer of the same type 
observed to be significandy increased in the selected rodent bioassay. These predicted 
relationships between effective dose and cancer ri.sk are then transformed into 
corresponding relation.ships between human exposure and cancer risk using an appropriate 
human PBPK mode! and a scenario for human exposure to CCI 4 in drinking water. Because 
almost all CKM models considered to date have assumed that premalignant cells may 
proliferate exponentially over time in vivo, this paper also examines the implications of 
using an alternative, perhaps more biologically reasonable, assumption'^ that such cells 
grow geometrically over time. Finally, a comparison is made between the "virtually safe 
dose" (corresponding to a lifetime cancer risk of 10~^) predicted using the CKM/PBPK 
approach illustrated here and that predicted using the more conservative approach typical in 
regulatory cancer-risk assessment. 

DATA ON CCU-INDUCED CARCINOGENICITY IN RODENTS 

Experimentally induced carcinogenic effects of CCI 4 observed in mice, rats and 
hamsters have led scientific and regulatory bodies to consider this compound a probable 
human carcinogen.Dose-response data on induced cell proliferation (see below) are 
only available for mice. Therefore, cancer-risk extrapolation for the present analysis will be 
based on cancer-bioassay data derived from studies on mice by the National Cancer 
Institute (NCI). In panicular, data will be used from the NC! studies^®"^^ of female B6C3FI 
mice weighing about 28 g that received 0, 1250, and 2500 mg/kg-d of CCI 4 , 5 days/wk for 
78 wk, by oral gavage using a corn-oil vehicle. Dosing began at 7 wk of age and surviving 
mice were observed up to 14 wk after dosing was stopped. Hepatocellular carcinomas were 
observed in one of the 80 pooled control female mice, 40 of the 40 low-dose mice, and 43 
of the 45 high-dose mice. 

Because 100% incidence was observed in the low-dose mice, these data cannot be 
used directly to derive a traditional estimate of hepatocarcinogenic potency for CCI 4 in 
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mice. Nevertheless, the data indicate that a dose of 35 mg CCI 4 per animal on the NCI 
dosing regimen yielded a cumulative risk of liver cancer greater than (40-1 )/40 or 97.5% by 
97 wk of age in these mice. The value of >97.5% is used as a lower bound on incidence rate 
in the analysis below. 


pbpk models 

The PBPK model of Paustenbach et al?^ was adapted to estimate uptake and 
metabolism of CCI 4 in mice and humans in the present study. Here, CCI 4 metabolism in 
rats was assumed to be described by a single saturable pathway with equal to 0.40 
mg/hr-kg or, equiv^entJy (by the body-weight scaling approach used here and described by 
Paustenbach et ai?^) 0.14 mg/hr for a 225-g rat (the value empirically determined^"^-^^ for 
rats cited, but then not used for simulation purposes^^). Corresponding scaled values 
used here for modeling 30-g mice and 70-kg humans are 0.034 mg/hr and 7.78 mg/hr, 
respectively. The present analysis m^es use of the partition coefficients for rats and 
humans used by Paustenbach et and assumes that mice have the same partition 
coefficients as the rats of the latter model. All other physiological parameter values for 25-g 
mice and 70-kg humans were taken to be the reference values reported by the U.S. 
Environment^ Protection Agency (U.S. EPA),^ except that the latter reported human 
^veolar ventilation rate was increased from 5.0 to a more realistic 5.6 liters/hr. Numerical 
implementation of the model was done as previously described.^^-^^ Other details and the 
results of PBPK analyses undertaken are described in the following sections. 

A study^^ of the absorption kinetics of four chlorinated organic compounds 
administered by a single gavage dose to rats in aqueous and com oil vehicles demonstrated 
that, over a 5-hr observation period, a com oil vehicle decreases the rate and extent of 
uptake of these compounds relative to that observed using an aqueous vehicle; the greatest 
decrease in uptake observed was for trichloroethylene, the least water-soluble compound 
tested. For the present an^ysis, oral absorption kinetics for CCI 4 were presumed to be 
similar to those observed for trichloroethylene administered in 10 mL/kg of com oil to 
400-g Wistar rats, which yielded a roughly constant, 90-min absorptive pulse of equivalent 
to a clearance of compound from the vehicle at a rate of about 0.11 mL vehicle/kg-min. The 
resulting total dose was about one 170th of that obtained using an aqueous vehicle (which, 
in turn, was assumed to be 100% of the administered dose). The vehicle-clearance rate of 
0.11 mL/kg-min used here was assumed to scale to the -0.26 power of body weight^^ and 
to apply only to single-dose exposure scenarios. The effect of using a com oil vehicle on 
absorption efficiency of chronically administered oral doses of similar compounds is 
currently unknown. 


DOSE-RESPONSE FOR INDUCED CELL PROLIFERATION 

Dose-response data on CCI 4 -induced hepatocellular proliferation in 30-g male CD-I 
mice were obtained from a recent study by Doolittle et The mice in this study were 
given single (or daily) oral doses of from 0 to 100 mg/kg CCLj in 10 mL/kg com oil, and 
measurements were made of the pciccnt of hepaiocytes in S-phase (%S) 48 hr after dosing 
and the concentration of selected enzymes in blood serum at various times after dosing. To 
estimate metabolized doses corresponding to the applied doses used in this study, one 170th 
of the ingested CCI 4 was assumed to be cleared from the com oil vehicle at the constant 
rate of 0.22 mLAg-min, yielding a 4S-min pulse injected into the liver compartment (see 
discussion above). The results of the PBPK analysis scaled to apply to the mice used by 
Doolittle et al. arc shown in Fig. la. 

Studies of rats expxised to CCI 4 by gavage have demonstrated that metabolism of this 
compound is the cause of its hepatotoxicity in rats.^*'^ One recent study further 
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1-hr meiaboliied dose (mg/kg liver) 


Fig. 1 , (a) Metabolized dose as a function of orally administered dose of CCI 4 in com 
oil to 30-g mice, based on the physiologically based pharmacokinetic model of Paustenbach 
et al. modified and parameterized as explained in the text. Solid and dashed curves 
represent 1-hr and 24-hr total metabolized doses, respectively. Solid points indicate the 
administered doses used in the study by Doolittle et al?^ referred to in the text, (b) Data 
(open circles) of Doolittle et al?^ on the percent of hepatocytes in S-phase 48 hr after a 
single gavage dose (or in some cases the fmst of daily doses) of CCI 4 in com oil to rats, as a 
function of the corresponding presumed effective ( 1 -hr total metabolized) dose derived 
from (a) above. Solid curve is Ae function given by Eq. (1) in the text, fit to the data by 
least-squares optimization. 


demonstrated that rats primarily metabolize CCI 4 to (Xij and found that the rate of CO 2 
production within 1 hr of dosing was the measure of CCI 4 metabolism most highly 
correlated (f?^= 0 . 88 ) with the observed degree of liver injury as measured by elevation in 
serum glutamate-pyruvate transaminase (SGPT) concentration.^^ It was therefore assumed, 
for the present analysis, that the total amount of CCI 4 metabolized in 1 hr is the most 
appropriate metric of effective dose, d, for the cytotoxic effects observed by Doolittle etal. 


The relationship between estimated values of d (from Fig. la) and the corresponding 
values of %S observed by Doolittle et al. is shown in Fig. lb. Also shown is the log-normal 
dose-response model 


%Si ^ 


logio(d) 

b 
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in which 4> is the cumulative normal distribution function, %Sq is the background level of 
(0.5% in the mice studied), and d is 1-hr metabolized dose. The model was fit to the 
(jata shown using least-squares optimization of the location and shape parameters a and b 
(yielding best-fit values of 0.23 and 0.13 mg/kg-hr, respectively). This dose-response 
function for CCIa-induced cell proliferation is a quasi-threshold model, based on the 
assumption that such proliferation represents regenerative cell growth following necrosis 
caused by acute toxicity of a sufficiently high dose of CCI 4 . The factor of 100%/3 in the 
model accounts for the fact that %S was measured by Doolittle et al. as a fraction of total 
hepatocytes undergoing replicative DNA synthesis at 48 hr post exposure, whereas not all 
hepatocytes induced into S-phase are involved in synthesis precisely at that time. Clearly, a 
synchronous induction of 100% of hepatocytes into S-phase would be lethal. Detailed 
experiments have shown that CCI 4 intoxication or partial hepatectomy induces in rats and 
mice a wave of regenerative hepatocellular proliferation corresponding to an S-phase pulse 
commencing between 24 and 48 hr ^ter administration and persisting between 
approximately 24 and 36 In fact, such proliferative response kinetics have been 

shown to be quite heterogeneous within rat liver, with hepatocytes in most liver zones 
reaching maximum %S values of between 20 and 35% at various times post injury.^’^'^* 
This is consistent with the estimated upper bound of 33.3% for %S used in Eq. (1), based 
on an average S-phase duration of about 8 hr (in rats)^® and an average duration of 
sustained induced peak S-phase activity of approximately 24 

The Doolittle et al. study included some investigation of how chronic (daily) dosing 
affects %S and serum enzyme levels in mice as compared to single dosing. Their results 
in^cate that %S is lower with daily as opposed to single CCI 4 dosing by a factor of 2 to 
3.^^^ For the present analysis, it was therefore assumed that the factor l(X)%/3 in Eq. (1) 
must be replaced by 100%/9 when this model is applied to chronic dosing situations. The 
latter assumption implies that no chronic effective dose, regardless of its magnitude, would 
produce a value of %S greater than 11.6%. 

An analysis of the Doolittle et al. data on CC] 4 -!nduced increases in %S and 24-hr 
SGPT concentration (Fig. 2) shows that the correlation between these two indices of 
hepatotoxicity may be sufficient to enable increased SGPT to be used as a quantitative 
predictor of increased %S according to the relation 


%Sd 


log 


10 


SGPTd'X 

SGPTo/ 


-h %So , 


( 2 ) 


in which the slope constant c is equal to 10.9%. Figs. 3a-3c summarize an investigation of 
this hypothesis using data on induced changes in SGPT levels in male 225-g rats riven 
single oral doses of CCI 4 in 2.5 mL/kg com oil from a study by Reynolds et a/.^^The 
PBPK analysis for this study (Fig. 3a) assumed that one 170th of the ingested CCI 4 was 
cleared from the com oil vehicle at a rate of 0.13 mL/kg-min, yielding a 19-min pulse 
injected into the liver compartment (see discussion above). The relation shown in Fig. 2 
was then used to derive predicted %S values corresponding to observed elevations in SGPT 
at the estimated effective doses of CCI 4 to rats for the Reynolds et al. study. Eq. ( 1 ) was 
then fit to these predicted %S values as a function of effective dose as described above, 
yielding the paranieter estimates a=0.93 and f>=0.12 rag/kg-hr corresponding to the dashed 
curve shown in Fig. 3b, where this fit is compared to the one pertaining to mice from Fig. 
lb. Fig. 3c shows that the dose-response models fit to the mouse data and the estimated rat 
data on elevated %S are approximately identical only if a from Eq. (1) for mice is nsmiTTwi 
to increase by a factor of 5 and if c from Eq. (2) is assumed for rats to be 5.45 (i.e.. half the 
value referred to above pertaining to mice). The first assumption suggests that the value of 
^Etw mice used in the PBPK analysis discussed above may be too low by a factor of 5, 
which is plausible in light of the fact that for mice in this analysis was simply 
assumed to be identical to an estimate obtained for rats (i.e., mice may have a higher 
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Fig. 2, Data (open circles) of Doolittle et al?^ on percent of hcpatocytes in S-phase 
(%S) 48 hr after gavage administration of various doses of CClj in com oil to 30-g mice, as 
a function of the fold increase (F) in SGPT concentration observed in these mice 24 hr after 
dosing. Solid line is the function given by Eq. (2) in the text, fit to the data by log-linear 
least-squares optimization. 


capacity to metabolize CCI 4 than would be expected based simply on allometric scaling 
from an experimentally bas^ estimate for rats). The second assumption suggests that the 
relation indicated by Eq. (2) and Fig. (2) may be species-specific. 


CKM MODELS APPLIED TO CANCER LIVER DATA FOR MICE 

The simple 2-stage CKM models used in the present analysis are described 
elsewhere and are reviewed in the Appendix. Specifically, the exponential and geometric 
CKM models used here to predict cumulative cancer risk, /?, as a function of age (minus a 
latency period here assumed to be 0 ) for a constant, lifetime exposure scenario, incorporate 
Eqs. (E-3) and (G-3) in the Appendix, where these arc replaced by Eqs. (E- 8 ) and (G- 8 ), 
respectively, for scenarios involving n distinct exposure periods (i.e., in cases involving a 
sequence of n sets of parameter values governing mutation and cell proliferation rates). 
Because the NCI bioassays for involved preexposure, exposure, and 

postexposure periods (see discussion above), Eqs. (E- 8 ) and (G-S) were used with n set 
equal to 3 to model hepatocellular cancer risk to the mice used. The cell-growth parameters 
of these models, g and a respectively (see Appendix), were fit to each of the single data 
points Rq= 0.0125 and R^f=€i.975 for the exposure scenario described above under the 
assumption that CCI 4 is a '’pure promoter," i.e., acts only to increase the proliferation rate of 
initiated cells. Here, d is the effective (1-hr total metabolized) dose to liver tissue for NCI 
bioassay mice, estimated to be 21.1 mg/kg liver using the PBPK model described above, 
scaled to 28-g mice under the assumption that 100 % of the chronically applied doses were 
absorbed. (Note that since 100% of these mice got liver cancer, Rj^>0.975 so the 
corresponding models actually describe a plausible lower bound on true risk. Fig. lb 
suggests that a similar response would have occurred had the effective dose been up to 
about 10 rinws less than it was.) To fit these models, it was also assumed that the NCI mice 
had L5 x hepatocytesfg liver tissue^ at risk throughout their lives of which roughly 
25% were diploid cells at risk for malignant transformation,^^ yielding a total of 6 x 10 
cells for 28-g mice with 5.5% liver tissue by weight, and that the two interstage mutation 
rales were each equal to 10 . The growth parameter values predicting the risks Rq and R^> 
iio<Sa) geometric and exponential CKM models respectively, produce 

the CKM functions shown in Fig. 4a. 
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Fig- 3. (a) Metabolized dose as a function of orally administered dose of CCI 4 in com 
oil to 225-g rats, based on the PBPK model of Paustenbach et al.^ as explained in text. 
Solid and dashed curves represent l-hr and 24-hr total metabolized doses, respectively; 
solid triangles indicate administered doses used in the Reynolds et al?^ study, (b) Open 
triangles represent the predicted % in S-phase 48 hr after single gavage doses of CCI 4 in 
com oil to rats extrapolated from Fig. 2 and from (a), using Eq. (1) from text; solid curve 
and associated data points redrawn from Fig. 1(b) for comparison, (c) The two curves in (b) 
converge if parameters arc changed as explained in text. 
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Incmased cancer risk 


Fig. 4. Exponential and geometric cell-kinetic multistage models fit to tumor 
response data^^^^ for female mice exposed to CCI 4 . (a) Solid and dashed curves are the 
functions 1 -exp(E) and l-exp(G), respectively, where E and G are given by Eqs. (E-3) and 
(G-3) in the Appendix, respectively, for which the constant XqitiiMx is set to 6 x 10“^. The 
functions are fit to the points /?o=0.0125 and R^O.975 at 97 wk, yielding the fitted 
exponential and geometric growth parameter values of [gQ=I.9I,g^6.54) and {aQ=fi.939 
and a^I2.7), respectively. The dotted line represents both models if the growth parameter 
is taken to be 0. (b) The fraction/correspxtnding to the growth parameter z = 

(where z is either g or a) used in the models from (a) to predict a given risk R, as a function 
of inert^ed risk for/fQ <R <Ri^. 


To proceed in the analysis, the background risk Rq was assumed to be independent 
such that increased risk is defined as 4 = {R—Rq)/{1—Rq), where Rq ^ R ^ and that 
premalignant (stagc- 1 ) clones are resistant to the cytotoxicity of necrogenic doses of CCU- 
but are stimulated to proliferate to the same extent that surrounding normal hepatocytes are 
during regenerative growth following dosing.^^ It was further assumed that a change in %S 
by a given fraction / of its observed range correspor Is to the same fractional change in net 
growth of premalignant cells, i.e., that 

f — ~ __ z’ — Zb ^3) 

^ ~ %Si-%So ~ Z4-Z0 ’ 
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where 0 </£ 1 and where zq < z < zj for z equal to g or a. Thus, by decreasing/from 1 to 
0 , the curves leading to in Fig. 4a are decreased until they converge to those leading to 
^Q. These assumptions and the information from Fig. 4a were then used to calculate/as a 
function of A for the exponential and geometric CKM models (Fig. 4b). Figure 4b shows 
that, based on the NCI mouse bioassay data, a reduction in/from 1 to between 5 and 20 x 
10 “^ is needed to ensure that/4=10“^. 


CKM MODELS APPLIED TO LIVER CANCER DATA FOR HUMANS 

The same approach taken above for using CKM models to describe liver cancer risk 
to the female NCI mice was applied to data on liver cancer incidence for human females. In 
this case, Rq was taken to be 0.0005, the cumulative risk of getting liver cancer by age 75 
for U.S. females during the period 1969-1971.'*^ For humans receiving the same effective 
dose of CCI 4 that the NCI mice received, the corresponding lower bound on increased risk 
was assumed to be approximately that applicable to the NCI mice. CKM models fit to 
each of these risk values are shown in Fig. 5a, and the corresponding relation of / to A 
shown in Fig. 5b demonstrates that, based on female liver cancer incidence data, a reduction 
in/from 1 to between 5 and 40 x 10^^ is needed to ensure that A=10~®. 


EXTRAPOLATION OF LIVER CANCER RISK TO HUMANS 

Figure 6 a shows effective dose as a function of / based on the data of Doolittle et 
air This relation was derived from Eq. (3) and the inverse of the relation specified in Eq. 
(1), appropriately modified to reflect a chronic dosing assumption. The analyses in the 
previous sections indicate that, at minimum, a reduction in /from 1 to between 5 and 
400 X 10 is needed to ensure that 4=10"”^, which, in light of Fig. 6 a, corresponds to an 
effective chronic CCI 4 dose rate of between 0.46 and 0.63 mg/kg liver per day. It was next 
as.sumed that this predicted range of "virtually safe" effective dose applies to humans as 
well as mice. To complete the analysts, a PBPK model for CC^-exposed 60-kg human 
females (with parameter values obtained as discussed above) was implemented using a 
chronic, daily scenario of ingestive exposure to CCI 4 , e.g., in drinking water. Exposures 
were presumed to occur once per day, each modeled to result in a constant 6 -min infusion 
of 100% of ingested CCI 4 into the liver compartment, and simulated exposures were 
repeated (for 7 simulated days) until dynamic equilibrium in effective (1-hr total meta¬ 
bolized) dose was approximated. The resulting relation of chronic applied to effective dose 
is shown in Fig. 6 b. This infonnation, combined with the previous analyses, leads to the 
prediction that an upper bound for the "virtually safe" dose rate of ingest^ CCI 4 for 60-kg 
humans is 3 to 5 mg/d, or 0.05 to 0.08 mg/kg-d, or 2 liters of water per day containing 
roughly 2 mg CCI 4 per liter. Recall that these represent upper bound figures because 100% 
of the low-dose NCI bioassay mice developed liver cancer. 


DISCUSSION 

The CKM-based approach to cancer risk assessment illustrated here for the case of 
CCI 4 differs fundamentally from more traditional risk assessment methods. Because the 
CKM model allows for a threshold-type dose response for agents considered to be "pure 
promoters," CKM-based risk estimates associated with exposure to such agents will 
typically be much lower than those based on "linear no-threshold" risk extrapolation 
methods. Uring the latter methods, for example, the U.S. EPA^ and the State of 
California*^’'*^ have, on the basis of animal cancer bioassay data such as the NCI study 
considered above, predicted a 10 "^ (upper bound) increased cancer risk to be associated 
with ingestion of 4 sym(180) 10“® and 5 x 10~^ mg/kg-d of CCI 4 , respectively, which are 
from 3 to greater than 4 orders of magnitude lower than the corresponding upper bound 
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Increased cancer risk 


Fig. 5. Exponential and geometric cell-kinetic multistage models fit to the liver 
cancer incidence rate'*^ {/?q= 0.0005) for 75 yr-old U.S. females during 1969-1971 and to a 
much larger risk (^,^.987) presumed to be associated with a high-level, lifetime exposure 
to CCI 4 . (a) Solid and dashed curves are the functions l-exp(£) and l-exp(C), respectively, 
where E and G are given by Eqs. (E-3) and (G-3) in the Appendix, respectively, for which 
the constant XQm|Afi is set to 1.8 x 10“^ (the highest value possible to ^low a Fit to the risk 
values mentioned. ITie fitted exponential and geometric growth parameter values obtained 
30^ (ao=00 and aj^=0.568}, respectively, (b) The fraction / 
corresponding to the growth parameter z = ZQ*f{z^Z(^ (where z is either g or a) used in the 
models from (a) to predict a given risk R, as a function of increased risk (/?-f?Q)/(l-/?^) for 


"virtually safe" dose rates estimated above using the CKM-based approach. Even if a 
reduction in the latter values by 1 to 2 orders of magnitude is made to account for the 1(X)% 
response in the NCI bioassay mice and for interindividual variability in humans, the CKM 
approach used here still predicts a much higher "safe" dose of CCI 4 than do the more 
tratlitional approaches. 

It must be emphasized, however, that large uncertainties remain in defining the most 
appropriate way to implement CKM-based cancer risk extrapolation, even assuming the 
CKM approach is comscL The most fundanaental source of uncertainty lies in the 
assumption made regarding the extent to which the observed carcinogenicity of a given 
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Fraction of original growth parameter ranga 



Fig. 6 . (a) Predicted cytotoxically effective dose to mice as the function of /defined 
by Eq. (3) in the text, based on the inverse of the function specified by Eq. (1) in the text (in 
which the constant, 3, is replaced by 9 to account for the effect of chronic dosing). Dashed 
lines represent the range in values of / corresponding to a 10“^ risk from Figs. 4(b) and 
5(b). (b) l-hr total metabolized dose at dynamic equilibrium as a function of daily oral dose 
of CCI 4 to 60-kg humans, based on the physiologically based pharmacokinetic model of 
Paustenbach et al.^^ modified and parameterized as explained in the text. Dashed lines 
correspond to those in (a). 


compound (say, in bioassay animals) is due to "promotion" (enhanced cell proliferation 
rates) alone, which reasonably can be presumed to have a threshold-type dose-response 
relationship, and how much is due to increased interstage transition rates (enhanced 
mutation rates), which is presumed by 2-stage CKM models to have a linear-quadratic 
dose-response relationship. Clearly, if a compound is assumed to have any degree of 
mutagenicity at all, this effect will ensure that the CKM dose-response relationship will be 
linear at sufficiently low doses.'** The analysis above assumed that CCI 4 is a "pure 
promoter," which, despite the very limited evidence for CCl 4 -induced genotoxicity in 
vivo, may not strictly be the case. Other sources of uncertainty in this analysis include the 
effects of com oil vehicle and chronic dosing on ingestive uptake of CCI 4 (sources not 
specific to CKM modeling), the selection of the most appropriate form of CKM model to 
use, and the quantitative relationship between changes in %S (or in SOFT levels) and 
changes in proliferation (birth minus death) rates in piemalignant cells. 
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Despite these uncertainties, the CKM approach incorporates the emerging facts 
regarding carcinogenesis to a far greater extent than do the more traditional risk 
extrapolation models. Thus, even when data available to implement a CKM-based risk 
extrapolation for a given compound are limited, consideration of this approach may provide 
heuristic advantages. For example, the present application to CCI4 suggests that levels of 
serum enzymes (such as SGPT) associated with cytotoxicity and/or induced cell 
proliferation be routinely collected from animals during cancer bioassays, and that 
application of linear or quasi-linear risk extrapolation methods to data obtained from such 
bioassays be ruled out for those dose groups exhibiting significantly elevated levels of such 
enzymes. 
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APPENDIX 

For cell-kinetic multistage (CKM) models in general,'^ the stage-; cell populations of 
a /t-stage process of carcinogenesis are represented by the stochastic processes 
j=l,...,k, over the period 0 £ f < T, where X^f0)=0; Xq(77 is a deterministic nonnegative (DN) 
("driving") function; and mfT), for ;=0,...,jl:-l and mj(TjT«l, are DN functions (stage- 
specific mutation rates) such that the probability of a unit gain by Xj+i(T) caused by a 
simultaneous unit loss from Xj(T) during the interval [T,T+dIj is approximately m;(T)dr. 
and the probability of more than one such event occurring in this interval is negligible. In 
semi-stochastic CKM models, wfTj, j—l,...JI:—l, are DN impulse-response (or weighting) 
functions used to represent (premalignant) cell proliferation subsequent to the occurrence of 
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each new stage-y clone arising by mutadon of a stage-f/’ - 1) cell, where = 1. Finally, 
let/f(f) = Prob(X^(r) SI), i.e., R{t) is the cumulative distribution of the waiting time for the 
occurrence by time t of the first increment in representing the first cancer cell. It has 
been shown for semistochastic CKM models in general that R{t) = for /?(f)«L and 

that R(t) “ l-€xp(£Y,j(f)) is a first-order approximation for all /?(?), where E denotes the 
expectation operator.!6 

Derived below are the approximate solutions for 2-stage (i.e, ;i:=2), semistochastic, 
exponential and geometric CKM models of carcinogenesis in the case that parameter values 
governing mutation rates and cell proliferation in intermediate-stage (here, in stage-1) cells 
change at times for i=l,...,n, where fQ^O, fj_i < r,, and t„=t. Here it is assumed that Xq{T) 
is equal to the constant Xq and that < T < /;) for J=0 and y=l are equal to the 

constants m,- and W,-, respectively. From these assumptions, it follows^*^ that for both 
exponential and geometric CKM modeh, 

” fti 

EXi(t„) / XorriiWiit, -T)dT , and (A-1) 

.= 1 

EX^itJ = -f f" M,EX,{T)dT . (A-2) 

In the following evaluations of Eqs. (A-1) and (A-2), the growth response function, wiin 
for all proliferating premalignani (i.e., stage-1) cells is taken to be a function of the growth 
parameter z, where z is designated g and a in exponential and geometric CKM models, 
respectively, and where z at time T equals the constant z^ for fj_j < Cj, i=Also, Pj 
will be used to denote the vector (m,vV/i,z,) of mutation- and growth-related parameter- 
values in operation over the period from /,_i to t=l, 

I. Exponential 2-Stage CKM Models 

For 2-stage exponential CKM models, wifF) = when «=1, so that Eqs. (A-1) and 
(A-2) yield 


£A'i(ti) % A'o —(e*-" -1) (E-1) 

9i 

= , and (E.2) 

EX,{i,) « -1) (E-3) 

= EX2{pi,t,) . (E-4) 


To consider cases in which «>!, we first define 


€r,.* 


( ~ ^i-1) for r < s 

i 0 for r > a 


(E-5) 


and note that wi(D = expfg^T+e i+i_„) during the period t^i to t„ for stage-1 clones arising 
at time T where < 7 < Upon substitution, Eqs. (A-1) and (A-2) yield 


EXi(tn) 


n -t,. 

St 


A’om; exp(gi(fi - r)-t- e.+i,„)dT 


(E-6) 


469 


PM3006451737 


Source: https://www.industrydocuments.ucsf.edu/docs/hkxj0001 



, and 


A:o—[ exp(Gi,„)-exp(e,+i,^) 


CE-7) 


» T + V [x,^^(exp(e.,.- 1 )) 

7^1 .<^n^ 9,gj 

• (exp(ejj -1)) expCei+,.j_i)] , (E-8) 

n. Geometric 2-Stage CKM Models 

For 2-stage geometric CKM models in the case that /!=1, it has been shown'^ that 
^l(D ~ {aT+iy (where in the text x=3 to model piemalignant clones assumed to grow 
spherically with radial growth proportional to time). For these geometric CKM models, 
^s. (A-1) and {A-2) yield 




F/2^2(11) ^ ATj 


{x -f- l)a 

= £‘.Y)(pi,ti) , and 

mjMi 


(t + 1 )(z -t- 2 )a^ 


(aiti + 1 )'^"- ((x + 2 )oi(i) 


+ 1 


(G-I) 

(G-2) 

(G-3) 

(G-4) 


To consider cases in which n>l, we first define 

7r, ~ I ^ + E,’=rO<(<t - <t-i) for r < .s 
1 1 for r > s 


(G-5) 


and note that ^((7) - during the period /^j to for stage-1 clones arising at 

timeT where rj_i < TSf,. Upon substitution, ^s. (A-1) and (A-2) now yield 


£A'’i(t„) ^ Y1 Xorrii (^a,(ri -T) + 7i-t-i,n) 


dT , 


(G-6) 


(G-7) 


EX,{i^) ^ Y. EX,{p„t,-U_,) + V [Xop-- 

•=i + ^){x+ 2)a,aj 


■ ( 7*+^ - - ■>•^+2 4 . ■v'^ + 2 'l 


(G-8) 


470 


* 



PM3006451738 


Source: https://www.industrydocuments.ucsf.edu/docs/hkxj0001 



